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Background 


¢ Google’s Earth Engine offers a “big data” approach to 
processing large volumes of NASA and other remote 
sensing products. https://earthengine.google.com/ 


¢ Interfaces include a Javascript or Python-based API, 
useful for accessing and processing over large periods 
of record for Landsat and MODIS observations. 


¢ Other data sets are frequently added, including weather and 
climate model data sets, etc. 


¢ Demonstrations here focus on exploratory efforts to 
perform land surface change detection related to 
severe weather, and other disaster events. 
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> Examples Classic “Hello World!” example... 








Online, shareable Interface for writing Javascript Output console for 
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print (aqua_ndwi_summer2010);_ 


Filter scripts... 1 // Create a 'collection', in this case Aqua MODIS Normalized Difference Water Index I) *ImageCollection MODIS/MYDO9GA_NDWI (91 e... Json 
2 var collection = ee.ImageCollection('MODIS/MYDO9GA NDWI'); type: ImageCollection 

> Private 3 // Filter that collection to the period of June 1, 2010 through August 31, 2010 id: MODIS/MYDO9GA NDWI 

> Shared (1) 4 var aqua_ndwi_summer2010 = collection. filterDate('2010-06-01','2010-08-31'); ST PEae 1449781080979000 
5 // Print the contents to the console a 

> Examples 6 >bands: List (1 element) 
7 


~features: List (91 elements) 
>O: Image MODIS/MYDO9GA_NDWI/MYDO9GA 005... 
>1: Image MODIS/MYDO9GA_NDWI/MYDO9GA_005 ... 
>2: Image MODIS/MYDO9GA_NDWI/MYDO9GA_005_... 
>3: Image MODIS/MYDO9GA_NDWI/MYDO9GA_005 ... 
>4: Image MODIS/MYDO9GA_NDWI/MYDO9GA 005 ... 


Sy aa p | st es mM p | om >5: Image MODIS/MYDO9GA_NDWI/MYDO9GA 005 ... 


+6: Image MODIS/MYDO9GA_NDWI/MYDO9GA_005_... 
1. Retrieve Aqua MODIS NDWI for all scenes in Summer (JJA) 2010. +i tenes upbin /unnesca wenn antesca 60s... 
2. Confirm retrieval by printing information to the console. gee carne tel ae a 


>11: Image MODIS/MYDO9GA NDWI/MYDO9GA 005... 
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// Create a 'collection', in this case Aqua MODIS Normalized Difference Water Index (NDWI) Use print(...) to write to this console. 
Addi ‘ var collection = ee. ImageCollection( ); 
D Filter that collection to the period of June 1, 2010 through August 31, 2010 
in // 2 . . f 
Get the same product 
n // Create a single image to display -- the maximum — each period. 
and period for 2011. var aqua_ndwi_summer2010_max = aqua_ndwi_summer2010.max(); 
var aqua_ndwi_summer2011_max = bask ial ganna t ect i 
Map.a r({aqua_ndwi_summer2010 max, ndwiViz, 
Compa as layers Map. aqua_ndwi_summer2011_max, ndwiViz, 
! r(diff, diffViz, )3 


var aqua_ndwi_summer2010 = collection. filterDate( ’ 3 
RCaueel ial the maxXIMmumM. // Difference these parameters 
// Display on a map. 
iterMap(-95,40,8); // Center display on the Nebraska/Iowa/Missouri borders 


On nuU PWN He 


// Create another period for Summer 2011. 
var aqua _ndwi_summer2011 = collection. filterDate( 3 
- var diff = aqua_ndwi_summer2011_max.subtract(aqua_ndwi_summer2010_max); 
Difference 2011-2010 
D | var ndwiViz = {min:0, max:1}; // Grayscale visualization for NDWI 
Disp ay var diffviz = {min:-0.2, max:0.2}; // Grayscale visualization for difference 
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1 // Create a 'collection', in this case Aqua MODIS Normalized Difference Water Index (NDWI) Use print(...) to write to this console. 
2 var collection = ee. ImageCollection( 'MODIS/MYDO9GA NDWI' 

3 // Filter that collection to the period of June 1, 2010 through August 31, 2010 
4 var aqua_ndwi_summer2010 = collection. filterDate('2010-06-01','2010-08-31'); 

5 // Create another period for Summer 2011. 

6 var aqua_ndwi_summer2011 = collection. filterDate('2011-06-01','2011-08-31'); 

7 // Create a single image to display -- the maximum for each period. 

8 var aqua_ndwi_summer2010_max = aqua_ndwi_summer2010.max 

9 var agqua_ndwi_summer2011_max = aqua_ndwi_summer2011.: 

10 // Difference these parameters 

11 var diff = aqua_ndwi_summer2011_max.subtract(aqua_ndwi_summer2010_max); 

12 


13 // Display on a map. 
14 var ndwiViz = {min:0, max:1}; // Grayscale visualization for NDWI 
15 var diffViz = {min: , Max: }; // Grayscale visualization for difference 





16 yer(aqua_ndwi_summer2010 max, ndwiViz, 'Aqua MODIS Max NDWI Summer 2010'); 
17 Map.addLayer(aqua_ndwi_summer2011 max, ndwiViz, ‘Aqua MODIS Max NDWI Summer 2011'); 
18 Map.addLayer(diff, diffViz, 'Difference NDWI Maximum'); 

15) c p ( ,40,8); // Center display on the Nebraska/Iowa/Missouri borders 
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Developer's Guide Python Insta lation ; ; . 
nstalling the client 
Introduction library 
Get Started The Earth Engine Python API is distributed as a Python package that is hosted on Github. The following instructions give Ubuntu Linux & Mac 
Earth Engine Code Editor an overview of installing the Google Earth Engine Python API. To use the Earth Engine Python API you'll need to instal! “ ——— 
Python Installation the client library and its dependencies on your computer and then set up authentication credentials. a 
Uninstalling the Library 
Image Uninstalling manually 
aneae ere Installing the client library hibetal 
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Image Visualization . 
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Ubuntu Linux & Mac OS X installation enn ee arenes 


Creating Images Coding in the Python API 


=e After the initial set up, the installation flows for Mac OS X and Ubuntu are nearly identical. 


Relational, Conditional and Boolean 


Operations 

Convolutions 1. Set up pip and Python Python packages are available 
Morphological Operations for the API aS well 

Gradients PIP is a package manager for Python. The following installation instructions assume that you are using it. 

Edge Detection 

Spectral Transformations hana tue Allows users to authenticate 
Texture Oo. 8 - - 

Object-based Methods Verify that you have Python 2.6 or 2.7: wa Interaction wide Ea rth 


Engine and write scripts, etc. to 
perform processing and other 
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ImageCollection python --version 


ImageCollection Information and 
Metadata 

sili allie cancion If needed, install 2.6 or 2.7 with apt-get . Then pip can be installed with: 
Filtering an ImageCollection 


Mapping over an ImageCollection sudo apt-get install python-pip Preferable in some Ways to 
Reducing an ImageCollection . f 
Compositing and Mosaicking USINE ce web client for 


Iterating over an ImageCollection 


interaction and coding. 





Mac OS X 


Geometry, Feature, FeatureCollection 
The installation instructions assume that you are using Mac OS X 10.9+, the Homebrew Mac OS package manager, and 


the pip Python package manager. Feel free to use a different package manager such as Fink or MacPorts if you prefer. 
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Applications: Land Surface Change 
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i} Docs Assets Disaster Case Studies/AprilTornadoesAnnualNDVI.js | Get Link | | Run ff} Reset ~Mxti foal Console ar. 
Filter scripts... 73 var ndviViz = {min: 0.0, max: 1.0, palette: ['000000', ‘O00CC00', '009900', '006600')}; Use print(...) to write to this console. 
74 var diffViz = {min:-0.2, max: 0.2, palette: ['CC3333','FF6666','FF9999','FFCCCC','CCFFCC','99FF99','6 
> Private 75 . -87.81 
~ Shared (1) 76 // Add the result to the map, assign the color palette, assign the name, tell it to display. 
~ SPoRT 77 Map.addLayer(pre_event_max.clip(polygon), ndviViz, ‘Pre-Event NDVI', true); 
A : 78 //Map.addLayer(post_2011_max, ndviViz, 'Max NDVI for May 2011', true); 
~ Disaster Case Studies # 79 Map.addLayer(diff_2011.clip(polygon), diffViz, ‘Difference 2011 and Pre-Event', true); 
hh 2011 August18HailMidwest E 80 //Map.addLayer(post_2012 max, ndviViz, 'Max NDVI for May 2012', true); 
& AprilTornadoesAnnualNDVi... oO 33 81 Map.addLayer(diff_2012.clip(polygon), diffViz, ‘Difference 2012 and Pre-Event', true); 
~ Sanity Check 82 //Map.addLayer(post_2013_ max, ndviViz, 'Max NDVI for May 2013', true); 
—gs. dipsiemmed : 83 Map.addLayer(diff_2013.clip(polygon), diffViz, 'Difference 2013 and Pre-Event', true); 
th MODISLandSurfaceTemperatures.js 84 Map.addLayer(diff 2014.clip(polygon), diffViz, ‘Difference 2014 and Pre-Event', true); 
& MonthlyMODIS.js 85 
& MonthlyMoDIS v2.js 86 // Center the map on the center of the polygon of interest 
~ Wildfi as 87 centerMap(center_lon, center_lat, 10); 
iatire 88 
hb daysSinceLastBurn.js 89 @ 
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} Docs Assets Disaster Case Studies/AprilTornadoesAnnualNDVI.js | Get Link | | Run ff Reset ~ Mu xtl ff... Console far. <| 
Filter scripts... 45 // Create a filter for the pre-event period. Use print(...) to write to this console. 
46 // Here, for the previous 31 days ending May 31 (i.e. month of May), years 2010-2010 (so just 2010). 
> Private 47 var pre = createPrevDaysFilter(5,31,31,2009,2010); -87.81 
~ Shared (1) 48 var post_2011 = createPrevDaysFilter(5,31,31,2011,2011); 
49 var post 2012 = createPrevDaysFilter(5,31,31,2012,2012); 
* SPoRT = 
. : 50 var post_2013 = createPrevDaysFilter(5,31,31,2013,2013); 
~ Disaster Case Studies /§ 51 var post_2014 = createPrevDaysFilter(5,31,31,2014,2014); 
ih 2011 August1 8HailMidwest i 52 var post_2015 = createPrevDaysFilter(5,31,31,2015,2015); 
& AprilTornadoesAnnualNDVI... Oo = 3 ee 
~ Sanity Check 54 // Obtain MODIS Aqua NDVI imagery within the polygon region 
ae a faceT . 55 var collection = ee.ImageCollection( 'MODIS/MYDO09GA_NDVI').filterBounds (polygon) ; 
anaSurtace Temperatures.js 56 
& MonthlyMODIS.js // Filter the collection to the time period of interest and obtain the max 
& MonthlyMODIS_v2.js var pre event_max = collection.filter(pre).max(); 
~ Wildfire // Now, the same thing for the post-event period (May 2011) 
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Filter scripts... 45 // Create a filter for the pre-event period. Use print(...) to write to this console. 
46 // Here, for the previous 31 days ending May 31 (i.e. month of May), years 2010-2010 (so just 2010). 
> Private 47 var pre = createPrevDaysFilter(5,31,31,2009,2010); -87.81 
~ Shared (1) 48 var post_2011 = createPrevDaysFilter(5,31,31,2011,2011); 
~ SPoRT 49 var post_2012 = createPrevDaysFilter(5,31,31,2012,2012); 
. . 50 var post_2013 = createPrevDaysFilter(5,31,31,2013,2013); 
~ Disaster Case Studies #51 var post_2014 = createPrevDaysFilter(5,31,31,2014,2014); 
fh 2011 August1 8HailMidwest i 52 var post_2015 = createPrevDaysFilter(5,31,31,2015,2015); 
& AprilTornadoesAnnualNDVI... Oo 7 5 ae ; 
~ Sanity Checks 54 // Obtain MODIS Aqua NDVI imagery within the polygon region 
Y : 55 var collection = ee.ImageCollection( 'MODIS/MYDO9GA NDVI').filterBounds(polygon); 
i MODISLandSurfaceTemperatures.js 56 
& MonthlyMODIS.js 57 // Filter the collection to the time period of interest and obtain the max 
& MonthlyMODIS_v2.js 2 var pre event_max = collection.filter(pre).max(); 
~ Wildfire ; 60 // Now, the same thing for the post-event period (May 2011) 
ib daysSinceLastBurn.js 61 a eee 
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| Scripts | Get Link Run Reset ~ m | Inspector | 
162 Click on the map to inspect the layers. 
163 

> Private 164 

> Shared (1) 165 

+ Examples 166 var mesh = ee. Image( '04525656613879846600-05414812628978701557'); qe NSS IVIESH 
167 var meshthresh = : 
168 
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Applications: Hail Damage 
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& Disaster Case Studies/201 x | Andrew | v 
€ C2 https://code.earthengine.google.com v = 
Google Earth Engine Search places and datasets... ai Sal Help ~ — andrew.molthan.work ~ 
Suce Docs Assets Porat gee Bi) IPA ae lea a Save icv nts Console Tasks 
Filter scripts... 218 var test = Chart.image.byRegion(ndvidiff, vectors, ee.Reducer.mean(), 250); Click on the map to inspect the layers. 
219 print(test); 
> Private 220 var test = Chart.image.byRegion(mesh, vectors, ee.Reducer.mean(), 250); 
> Shared (1) 221 + print(test); 
> Examples 222 


223 Map.addLayer(vectors, {}, ‘hail vectors’); 
224 Map.addLayer(hailmask.clip(polygon),{},'hail mask’); 







BYO-Data: 

Imported data sets can be used alongside other Earth Engine 
holdings to provide additional information for algorithms, or 
Earth Engine techniques can be applied to input data sets. 
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Based upon technique explored by Bell (2015), discussed on Thursday, 2:30, rm 252 


Strengths / Weaknesses 


¢ Strengths 
¢ Earth Engine data holdings are tremendous, and the streamlined 
Javascript/Python API is capable of processing large volumes of 
data very quickly. 
¢ Relatively simple, big data platform for exploring various remote 
sensing products and capabilities. 


¢ Numerous and increasing number of nigh volume data sets 
available that could be difficult to download and process locally. 


¢ Weaknesses 
¢ Primarily focused on raster types of data — large gridded or swath 
products of comparable resolution. 
¢ Most products still require some understanding of the raw data — 
many appear to be a literal depositing of the source data with no 
additional post-processing. 
¢ Ex: Temperatures in scaled units, not K 


¢ Documentation could benefit from a few use examples, and some 
experiences where documentation or errors were cryptic. 


Questions? 


¢ Check out some related talks! 


¢ Near-Real Time Severe Weather Damage Identification Algorithm 
for Vegetation: Development and Early Results 


¢ J. Bell, University of Alabama in Huntsville 
¢ Thursday, 3:30 pm, 252 


¢ andrew.molthan@nasa.gov 


